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Introduction 


V/Hiv-l Elni .a element methods can be considered as con- j 
/enient methods for solving partial differential equa- | 
jtions. In the case of solid mechanics the most convnon ' 
jf°naulation is to use displacements as field variables ! 
JM ! 1n the d1f Eerential equations. The so-called conven- j 
j t10nal finite element formulation is an approximate ! 
§5§ (solution based on assumed displacements and the corres-' 
£Ii | pondins varia tional principle is the principle of tne * 
Wi f Stat1onap y Potential energy. Such method can thuSj I 
WM jbe called a Primal Finite Element Method. For solid ! 
SS mechaMcs Problems a dual formulation can be made using 1 
ggf .stress functions as field variables and the correspon- ; 

HI ! d<n9 Comp1enientar y en ’ ep sy principle. Applications of ‘ 
j 806 * 1 PHincIple for the finite element method have been 
*made, but' are not popular. j 

HJ At a rather ea:-1y age of finite element develop. 1 
||g ment, short-comings of the assumed displacement method j 
were discovered. The main one Is the difficulty In 1 
1§! constructing shape functions for thin plates and shells- 
Iff which require C continuity-conditions. Other obs^rva- 1 ; 
H| t,ons are that the so-called compatible elements are 1 
t0 ° rl9<d and * * ndeed » for some lilting cases, such as : 
S wly Incompressible materials or thin plates an<T : 

“M 5he11$ Wh1ch are form «lated by talcing transverse shear-^ 

! P n9 stra ' ns Into account, the corresponding finite ele- 
ment model may be completely locked. Another area that 
icalls for Improvement is to avoid the loss in accuracy 
«. 2i - - when the ev3luat1ofl of strains and hence stresses are 
i§; acc emplished by differentiation of the assumed 
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In two new formulations of Hybrid/Mixed finite 
element, methods respectively by the Hellinqer.Re1s«ner 
principle and the Hu-Washizu principle, the stre« equl,' 

11bp1um editions are brought in as conditions of con- j 

stra1nt through the introduction of additional internal' 
p| splaeement parameters. The new approaches are more j 
flex1ble and hav e better computing efficiencies. A j 

•^procedure for the choice of assumed stress tenns for ♦ 

g||j3-0 solids is suggested. Example solutions are- given? 
^||| |for plates and shells using the present formulations <‘ 
mWavHv jAbd the idea of semiloof elements. I 


displacements. [ 

Many alternative finite element methods that have [ 
been developed are bared on the introduction of addi- [ 
. tiona) field variables. Within the element more than ! 
one field variable can be introduced such as methods 
! derived by the Hellinger-Relssner principle for which 
(stresses and displacements are used and methods derived- 
Iby the Hu-Washizu principle for which stresses, strains!’ 
land displacements are used. Such models are named f 
jmixed models. .When additional variables are introduced- 
jas Lagrange multipliers for the maintaining continuity 5 
conditions along the mterelemen: bouncary it is la- I- 
belled as a hybrid model. These two models are not 
mutually exclusive and in many cases the same finite j 

element method can be derived using different varia- [ 

ticnal formulations. It is, thus, decided tha\ a dis- f 
cussion of non-primal finite element metnod be Vumced ! 
'together under the name Hybrid/Mixed finite element 
| method. . j 

j The assumed stress hybrid method which makes use \ 
t of boundar y displacements as additional variables, and f 
(hence leads to matrix displacement methods in the fi- ! 
jnite element formulation, can overcome the shortcomings!' 
j°f the assumed displacement methods that were mentioned-' 
.above. However it also needs further Improvements. * 

; The main ones are (1) the lack of a guideline for the i 
‘choice of the assumed stress terms, ( 2 ) tne excessive L : 
.computing time needed in the inversion of a flexibility' 
j matrix during the process of constructing the element' f 
I stiffness matrix and (3) the difficulty in the choice f 
jof equilibrating stresses in the case of complex prob- «' 

| lems such as shells. | 

j . The objectives of the present lecture are: E 

i(l) To present two recently developed Hybrid/Mixed •! 

(finite element formulations based, respectively on the \ 

; Hell inger-Reissner principle and the Hu-Washizu prlnci--’ 
.Pie for which stress equilibrium conditions need not be*' 
(Imposed Initially but* are brought In as conditions of \ 
constraint, instead, by Including additional Internal 
•displacement parameters. These new approaches make the 1 
.formulation more flexible and also improve the comput- : 
Ing efficiency. \ 

(2) to present a systematic procedure for the choice 
of the ^ + tha* will 
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v^ ; ’prevent the appearnace of kinematic deformation modes, 
i;:!;;’: Three dimensional elements are used as examples. 

'( 3 ) t 0 present some results on plates and shells using 
iresent approaches and the Idea of semILoof elements.!*? 


, 2 Hybrid/Mixed Formulation by Hellinger-Relssner 
Principle 

Element stiffness matrix can be constructed using 

'different variational mothods. The Initial step In com- 

(non Is to assign nodal displacements q and correspon- 
mi;' : ding boundary displacement u which maintain compatlblll- 
V/.V.V. ty with neighboring elements. Then ty using any varla- 
tlonal principle In solid mechanics a solution for such 
' a problem with prescribed boundary conditions can be 
obtained. Such solutions may be in terms of displace* 
jments, stresses and/or strains. The element strain ; 
’energy U can thus be expressed in terms of the nodal 
^displacements q, and the element stiffness matrix k can 
fa* obtained. * ; 

iiijjii. | * In using the He11inger»Reissner principle to . 
;0}} .solve the above boundary value problem, the following 
functional should be maintained stationary, 


"R 


/ c- 

v 


r-0 T S 9 * o'CD u)] dV 


( 1 ). 




stresses o then the resulting element is Identical to 
'the compatible element, as was pointed out by Fraeljs 
•de Veubeke.^ If o satisfy the homogeneous equilibrium, 
^equations 

t 

j ♦ D T o - ° (3) 

* 

jThen the second term In Eq. (1) can be changed to the 
boundary Integral, resulting in 


*p • /[- 4 o T S o] dV + / T'u. ds 

n y £«►**«• ^y - •* 


rTn 




(4) 


1 


This Is the negative of the element complementary ener- 
gy and Is the basis for the Initial development of the 
:assumed : tress hybrid element,^ which may be consi- 
dered as the result of a hybrid variational principle. ^ 

:Th1s also indicates that if the assumed stresses o In 
* ** 

Eq. (1) satisfy Eq. (3) and u a u the resulting element 
stiffness matrix by a mixed variational principle Is 
Identical to that bv a hybrid variational principle. 


•2* 


- / T (u - u ) dS 
j 3 V * - “ 

• : 

‘where Sis the elastic compliance matrix and 

i I • i 2" (2)' 

a 

jrelates the boundary tractions and the stresses. Here 
the element displacements are assumed to be. In general, 
not compatible with the prescribed boundary dlsplace- 
;ment u. If the element displacements u are compatible 
•with u, and there Is no constraint for the assumed 


With stresses expressed In terms of stress parame- 
ters e, 


j'N 


and the element displacements u or boundary displace- 
ments u interpolated In terms of nodal displacements 3 , 
l.e. 



U » N 3 


(6) 


U, U 


(7) 

then 

ff c * * ’’r ’ I S 

- B^G 3 

(8) 

where 

H . / p T S P dV 

* y * • • 

• 

(9 j 

and 

G ■ / P T N dV by 

V * * 

*R 

(10) 


■ / pV^Nds try 
av 

7r c 

dlj 


k « gV 1 G 


( 5 ) 


Then 6it » 0 yields 6 In terms of q. The resulting ele- 
ment stiffness matrix obtained from the strain energy 
expression Is ! 


( 12 ) 


In a new formulation, the stress equilibrium con- 
ditions Eq. (3) are not satisfied but the element dis- 
placements u are expressed in terms of nodal displace- 
ments q and additional internal displacement parameters 
X , 1 . € * * 

U = u q +u x - Nq + MX (13) 

Here Uq are compatible with the boundary displace- 
ments u but u x are introduced so that u are no longer 
compatible. By Integrating the term with u x and by 
recooni 2 lng that 

# 

u. » u - u on 3V 

-X - - , 

one obtains from Eq. (1) ; 

ff R ■ {t4 2 + 2 T( 5 Hq) * <3 T 2> T <*V 04) 

If are bubble functions, l.e. u x along the 
element boundary are zero, then u - u^« 0 and the above 
equation still holds. It is seen that the last term In 
Eq.(2)conta1nS the stress equilibrium condition with u x 
'as Lagrange multipliers. Thus, In finite element for- 
mulation the stresses heed not satisfy the equilibrium 
condition Initially while the Introduction of the in- 
ternal displacements u x will enforce this condition. 

With 0 , and u x represented by Eqs. (5) and (13), 

*R * “ 7 B T H 6 + B T 6 3 - g T R X (15) 

where H and G are given In Eqs. (9) and (10) and 

R ■ / (D T P) T M dV (16) 

The stationary condition of ir R then yields 
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6 « H*' (6 q * R X) 


Rig « 0 

VtVAV.VtV/tt*. * * 


HXi::::; Eliminating X and substituting 8 Into the strain energy 
K&wwyltxp re s s 1 0 n one obtains 

jtt 

88§ k « gV’g - gVr (rVr)’ 1 rVg (19) 

m m m m % m m m m m mm m 

Equation (18) Is, of course, the equilibrium constraint 
'Eiiiiiii * or the stress parameters 8. When problems are set up 
»it"f;:::;;Vjln terms of Cartesian coordinates, such constraint equa- 
tlons can be Identified directly, for thin plates and 


shells for which compatible shape functions are diffl- 
cult to construct, * c of Eq. (4) Is used by the exist* 
T££&"jing hybrid stress method. Tong has suggested the use 
L;:HL of uncoupled stresses while adding a constraint term 
L to Eq. (4). 6 The corresponding variational func* 


tional may be written as 


I *d * /C* i c^S o - (D^o)^ u,] dV + / T T u dS 

1 V 3V (20) 
'and the resulting element stiffness matrix Is the same 
as Eq, (19). The motivation by Tong is that when the 

2 matrix in Eq. (5) are not coupled, the inversion of 
the flexibility matrix H can be simnlified ^*®. For 
‘example, in the case of 3*D solid, with 5 * {o x o y o ? 

•°xy c vz C Kz' i ’ ^ the exSans '' on f° r stress 
■components are of iaentical form then the H-matrix 
•Eq. (9) will take the form 


where S Is the elastic compliance matrix and 


t • 1 Ei Ei a » 

The inversion of H then becomes simply 


,41 «*.* C % 


where C » S*^ and £ ■ . (24) 

Although In using uncoupled stresses the size of 
H Is Increased and the inversion of another matrix must 

* <w 

be performed when Eq. (19) Is used to derive the 
element stiffness matrix, there could be considerable 
savings In computation when the orders of <|i are much 
■smaller than that of the H matrix of a coupled stress 

* I t h •**♦-«* ♦Kef <r\+ 


are coupled In the stress-strain relation It Is only ? 
necessary to use the same stress expansion for the t 
normal components. ;• 

The present formulation by using additional 
Internal displacement parameters Is a more general way : 
to Introduce the stress equilibrium conditions. For 
exampla, the stress terms may be expanded In terms of 
natural coordinates such as Isoparametric coordinates 
for qualri lateral and hexahedral elements and area and • 
volume coordinates for triangles and tetrahedrons. 

Also, a 'imlted number of additional terms may be used. 
In these cases the stress equilibrium equations are 
, satisfied approximately but the formulation of the 
problem Is simplified. It should be remarked that In 
the formulation by the Helllnger-Reissner principle ; 
the equilibrium conditions need not be satisfied for 
the assumed stresses. . , j 

3 Hybrid/Mixed Formulation by i 

' Hu-Washlru Principle . j 

In using the Hu-Washlzu principle to solve the 
above boundary value problem the following functional 
should be maintained stationary, 

it hw 8 / [4 c T C c • c T e + o T (0 u)] dV 

- ; T T (u - u) dS (25) 

3V ' “ 

In 1973,' Wolf^ used the Hu-Washizu principle for 
finite element formulations, but the resulting method 
has not been used widely. l.K. Tang and his 
colleagues^® at the Dalian Institute of Technology 
have initiated the so-called quasi-conforming elements 
or string-net elements for the construction of element 
stiffness matrices and have discussed the connection 
of their method with the Hu-Washlzu principle^. An : 
Interpretation of their method is that in applying : 
the Hu-Washizu principle each stress component and j 
Corresponding strain component are expressed In terms j 
of the same polynomial expansion, and the rasultinq ; 
matrix to be Inverted then becomes one with only i 
diagonal submatrices. This approach can again be 
combined with the use of additional internal displace- 
ment parameters to enforce the stress equilibrium 
■ conditions. Aqain, using the 3-D problem as an 
;* example, In the finite element implementation, let 

c • P a (26) 


with Identical P for stresses and strains but the 

m 

expansions for different components may be entirely 
Independent* With element displacements again 
evresned bv Eq, (13), equation (25) then becomes 


t 4 


) 


tv* ** f 


oiidMo. tkx is 

Of focm QUALITY 

r Hw * I 2 T j! * * s T § 2 + S T g 3 * c t R X ( 20 ) 

(29) 

H< ' ( 30 ) 

-»G 


where J ® / P T c P dv 

*> y *» *» *» 

H • / P T P dV • | H 

f * y *» •* 


1. 


and G and R are given In Eqs, (Id) and (16), 

• Th0 stationary condition of r HW with respect to 
8 a and X then leads to three sets of equations, which 
■can be used to obtain a In terms of g, The resulting 
.element stiffness matrix Is 


G'W G - G T K R (R T w R)* ] r t w G 


(31) 


l k « G T W 

» «* 

I 

- where « . fl-'j J*1 , 32) 

As shown In Eq. (30) H matrix consists of only diagonal 
submatrices. Thus, Its inversion always involves only 
the Inversion of H,. even if the material property is 
anisotropic. ‘ 

For general quadrilateral and hexahedral elements 
It Is more convenient to express strains, stresses and 
displacements in the Isoparametric coordinates. In 
that case, the Jacobian !J| appears in the area or 
the volume integrals, here, however., the strain terms 
can be expressed as 

• J l £ * t 3 (33) 

and for properly chosen P, the resulting submatrices 
may all jecome diagonal matrices and R will also be 
a simple diagonal matrix. 


4 Choice of Assumed Stress Terms 
for Hexahedral Elements 

The guidelines for the assumed sfesses in the* 
hybrid/mixed finite element formulations are to include 
aoequate terms to suppress all possible kinematic 
t Reformation modes and to avoid excessive stress terms 
' In order to prevent the forming of an overly rigid 
element. For hexahedral elements of general distorted 
geometry kinematic deformation modes can, In general, . 
be -suppressed when the number of stress terms Is equal 
or larger than the total number of degrees of freedom 
minus six. However, for regular brick-shaped elements 
The possible displacement modes are of simple shapes, 
not all assumed stress terms are useful for suppress* 
ing the kinematic modes. Thus, their selections become 
a very Important step In the finite element formula- 
tion. 


-4- 


Now consider a 20-node brick element with the x, 
y and z axes coincide with the axes of symmetry ns 
shown in Figure 1. The element displacements can be 



Figure 1 20-Node Brick Element 
expressed in terms of the nodal displacements 


rff (1 ' 

1 to 20) 

as 
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Ki(x>y>z>) - 

v i ■ 

w f j 


u i V 1 


(34) 


One can conclude that the shape functions N } (x,y,z) 
can be represented by 20 independent and simple 
polynomial terms as follows: 

Constant, x, y, z, x^, y*, z?, xy, yZ, X 2 , x^y, 

^ ^ # x i, kz « rX/ 2 , .x^yz, x y^z and xyz^ 

For a regular brick element, then, all possible element 1 
displacements can be represented by 6C independent j. 
displacement nodes. Here, the six rigid body displace- 1 
ment modes correspond to zero strains, hence, can not l 
be suppressed by any stresses and one only needs to 
consider the remaining 54 displacements. 

It is seen that each displacement mode may 
correspond up to three non-zero strains. Let e be 
one of these strains then, to suppress a particular 
mode it Is necessary that a stress term o exists for 
the stress component corresponding to e such that the 
integral / y ocd V does not vanish. Each stress term, 
however, can be used to suppress only one Individual 
displacement mode. An ideal situation, thus, would be 
to choose 54 assumed stress terms each of which Is 
associated with one individual displacement mode. 

One of the possible choices is shown In Table 1. 
Here, In the various block spaces the polynomial 
expansions for the six stress components and the 
corresponding displacements to be suppressed are 
indicated. There are 54 Independent e-terms. In many 
other blocks stress terms are Introduced to maintain 
the equilibrium conditions. They are not used to 
suppress any of the 54 displacement modes. 

It should be noted that some displacement terms 
may each correspond to three strains and some strain 
terms may be arrived from more than one displacement 

mode. The net remit •(» *,.» e* 


t 

I 
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f-VV: 
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: Table 1 

j Stress Terms and Corresponding'Displacement Modes 
i to be Suppressed for 20-Mode Brick Element. 
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C 1 

u*x 
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V*y 

h 

w*z 

'4 

u B y,v»y 

*6 

vz.vpy 

-6 

^Z.W’X 
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**14 

*®21 

e 9 

v*xy 

*11 

W n X2 


*16 
w «xy 

*19 2 

w®x 

y 

h 

u*xy 

* 6 13 

* 6 1P 

*12 

vryr 

% 

u*y 

6 17 2 

w*y 

e 20 

u°yz 
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*8 

U»X2 

*10 

v°y2 

**17 

*h? 

*15 

v»xz 
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V®2 
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U®2 

xy 

e 22 

u»x 2 y 

*27 

v»xy 2 

*32 
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t 4l 

w*x 2 y 

yz 

6 23 

u«xy 2 

*28 
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h, 
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*37 

2 
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X2 

*24 

u=x 2 z 

*29 
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*34 

w*xz 2 

*38 
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v*x 2 y 

*35 
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w »x z 
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y 2 

*25 2 
u*xy 


*35 
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•v 2 
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displacement modes may be suppressed by alternative 
stress terms and there exist many other combinations of 
B4 stress terms that can suppress all possible 
kinematic modes. 

Another observation Is that the stress expansion 

Is complete only up to linear terms, hence, the element 

properties will be dependent on the reference axes fpr 

the assumed stresses. Spllker ot al . 12 has shown that 

a hybrid stress element may behave excellently under 

one coordinate system but will contain a kinematic 

deformation mode under another reference system. For 

a regular brick element, the reference axes should 

obviously be the axes of symmetry. For elements of 

general geometric shape it Is still required to choose 

a local reference coordinate system In order to achieve 

an optimal element behavior. . , 

It is also noted that there are alterative 

choices of the stress terms for maintaining the stress 

ecui librium conditions. In Table 1 complete symmetric 

conditions are not maintained In this particular stress 

selection. If such conditions are to be maintained 

ard all the stress equilibrium’ equations are to be 

satisfied the number of independent stress terms will 

necessarily exceed 54. Alternatively, some stress 

eauilibrium conditions can be relaxed in order to 

maintain complete syrretric conditions. It can be 

seen that in Table 1 if we exclude all the underlined 

terms which are Introduced to maintain the equilibrium 

conditions of all cubic stress terms, the symmetry 

condition is satisfied. An example solution of a i 

sphere with a hole at its center subjected to varying ! 

temperatures along the radial direction shows, in factj 

that by using fewer stress terms in such a way the J 

results are better than those obtained by using all f. 

terms in Table 1, and even better than that by the 

conventional displacement method. It has been found, 

however, that the equilibrium conditions for the lower $ 

order stress terms should be maintained in order to 

achieve an optimal element behavior. Also excellent j’:; 

results can be obtained by using an element with 57 

£*terms obtained by replacing the e 52 , gj 3 and 

tenns in Table 1 by o- s^xy 2 + e S 3 xz 2 , a • e^x 2 y € 

+ ^ 55 ^ * °z * 2 + ® 57 x * 2 without any shear 

stress terms for maintaining equilibrium conditions, fe 

Details of these solutions are given in Ref. 13. bs> 

For brick elements with 8 , 12 and 16 nodes the 0;; 

independent displacement modes can all be identified 

from the 60 modes of the 20-node element. Thus* the 

corresponding stress terms can also be identified jv.vv 

immediately from Table 1 , For example, the indepen- \rZ‘. 

dent displacements for the 6 -node brick element are i;;c ; 

r .**.*• 
t;v; 

f:V*; 
l *.\v 


u ♦ v » w * constant, x*y, 2 ,xy,yr,x 2 ,xy 2 
and the corresponding assumed stresses are 
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° x r e 1 * e 2 y + 6 3 * + ;y* 

°y " e 5 ♦ figX ♦ B ? * + BgX2 

°* ’ p 9 + 6 10 x + *11* + <V* 


°xy * °J3 * *14* 


‘yt " P 15 + ^g x 
°x* * *17 + *18* 


(36) 


In this case the three quadratic stress terms are 
, necessary in order to suppress the displacement modes 
u, v, w ■ xyz. It Is also noted that even the linear 
terms arc not complete here. In fact, If the linear 
terms are complete with altogether 24 0's, the result- 
ing element will be almost as rigid as an assumed- 
displacement element based on trlllnear shape functions. 

In using the present methods by r R for Isotropic 
materials, the stress expansion of the three normal 
stress components must be Identical while the Shear- 
ing stress components may bo kept the same as those 
In Table 1. For a formulation similar to the £7-0 
! element described earlier, the size of the H-natrlx 

I «* 

; is 75 x 75 and the evaluation of an element stiffness 
! matrix will require the Inversion of one 18 x 18 
j (SV 1 matrix, one 17 x 17 c-matrlx and three 8x8 
j matrices instead of the Inversion of a 57 x 57 flexi- 
bility matrix. When the formulation is by and the 
stress terms are the same as that by r R , the result 
will be identical to that by r R . However, in this 
case the stress expansions for different normal 
Components heed not be the same and one may use fewer 
stress terms than that by n D . 

K J 

* I 

5 Serr.lLoof Elements for Plates and Shells 

. i 

For thin plates and shells* Irons has pointed . 

; out the advantage of using semiLoof elements^* for 
| which the derivatives of lateral displacement w are 
| not maintained compatible at the corners of the 
! elements while normal derivatives w at some points 
j along the interelement boundary are kept the same* 
j 8 y using the assumed displacement approach the derive* 

: tlon of a semiLoof element is very complicated* 
j However* it is rather straightforward using the assumed 
| stress hybrid method 12 . Within a semiLoof element the 


derivatives of boundary displacement w are not contin- 
uous at the corner nodes but the interelement compat- 
ibility conditions are completely satisfied. 

Triangular plate elements with 12 D.O.F. as shown 
in Figure 2, have been derived using 9, 13 and 17 
I Independent moment terms. The Loof-nodes are located 
j at 1/3 points. The 9 moment terms are the complete 
I linear terms, hence, resulting element by the hybrid 
, method is identical to the equilibrium element by 
^ Prapi 1* it 


*6- 




Fig. 2 Center Deflection of Clamped Square Plate Under 
, a Concentrated Load 

t 

♦ 

.as shown by the solutions (Fig. 2) of .the center deflec- 
tion of a clamped sauare plate under a concentrated 
load. The use of four additional moment terms improves 
the accuracy of the solution. The addition of eight 
moment terms, however, makes the element too rigid. 

'Two 16-dof semiLoof elements of rectangular shape have 
also been formulated using 13- and 17-6 terms. These 
elements are necessarily more flexible than the corres- 
ponding triangular elements as shown by the results • 
given in Figure 2. I 

For shell elements the equilibrium equations will 
always couple the membrane stresses and the moment 
stresses. However, for Isotropic materials, they are 
not coupled in the stress-strain relations. Thus, in 
this case, by using the .new approach by the H-matrlx 
can always be partitioned into four sub-matrices along 
the diagonal. In the new approach by * HW , of course, 
different components In the membrane and moment stresses 
are always decoupled. Another added advantage of the 
new .formulations Is that the resulting H matrices by 
both variational principles do not contain the geometry 
of the shell, making the implementation simplified. 

Both triangular and quadrilateral shallow shell 
elements have been constructed for semiLoof boundary 
nodes by the present methods based on both the 
Hell Inger-ReiSsner principle and the Hu-Washi 2 u princi- 
ple . The assumed stress terms which for the two 


r: 

l 
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'v:::Zv. Fig. 3 Stress Distribution of a Pinched Cylinder by 
Mfi-'l. Semi Loo, f Elements • : 

•L".' i 
€ lenient* cbnsist of 37 independent £'$ initially* and 
AUIUI;-::, 12 constraint equatic for £ are then introduced 
reducing the number of independent B's to 25* For 
the 32-dof quadrilateral shell element 55 independent 
6 1 s were used initially and 22 constraint equations 
Safes 2. for 6 are then 'introduced reducing the number nf 

Independent 8's to 33. I 

f§£ A problem of the pinched cylindrical shell with , 
freely supported ends was analyzed using 4x4 mesh. 
p-Mk It was found that the displacements and stresses 
obtained by these two variational methods differ so 
WIM; HtH* the resulting plots of the membrane and 
frfrfjiii moment distributions (Fig. 3) can not be distinguished 
ij: Z!:;ik from one to the other. An analyses of the computlhg 
efforts for these two approaches Indicates that the 
JI’MM formulation by’ the Hu-Washlzu principle has an edge 
over . that by the Helllnger-Relssner principle. 

methods are simpler to formulate than the original 
hyorld stress method for which the membrane and moment 
••V'rrV;' stress terms are coupled by the equilibrium equations. 

' 6 Conclusions I 

• i 

* » 
h;;£± New applications of the Helllnger-Relssner j 


principle and the Hu-Washi^u principle to finite 
element formulations have led to a more flexible and Pj~ 
most efficient finite element method by the assumed 
stresses. The possibility of using natural coordinates 
and the relaxation of some equilibrium conditions are 
the flexibility of the new hybrid/mixed stress method. 

The decoupling of the flexibility matrix H leads to a * 
saving In computing time. A systematic procedure for 
the choice of stress terms In certain hybrld/mixed 
| finite elements has been established based on the 
requirement for suppressing the kinematic deformation [ 
modes. Numerical results indicate the advantage of | 

; of maintaining the equilibrium of the lower order j..".." • 

* • f 

J stress terms while relaxing that o’ the higher order 

t terms. Hybrid/Mixed formulations are suitable to 

J plates and shells using scmiLoof elements. j 

i ' i 
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